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Abstract 

The difficulties arising in the investigation of finite-size scaling in d-dimensional 0(n) systems with 
strong anisotropy and/or long-range interaction, decaying with the interparticle distance r as r~ d ~ a (0 < 
a < 2), are discussed. Some integral representations aiming at the simplification of the investigations 
are presented for the classical and quantum lattice sums that take place in the theory. Special attention is 
paid to a more general form allowing to treat both cases on an equal footing and in addition cases with 
strong anisotropic interactions and different geometries. The analysis is simplified further by expressing 
this general form in terms of a generalization of the Mittag-Leffler special functions. This turned out to be 
very useful for the extraction of asymptotic finite-size behaviours of the thermodynamic functions. 
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I. INTRODUCTION 



The standard finite-size scaling theory (FSS) is usually formulated in terms of only one ref- 
erence length - the bulk correlation length £. For a system with finite linear size L, the main 
statements of the theory are: 

(i) The only relevant variable in terms of which the properties of the finite system depend 
in the neighborhood of the bulk critical parameter (the temperature in classical systems and the 
corresponding quantum parameter in quantum systems) driving the phase transition is L/£. 

(ii) The rounding of the thermodynamic function exhibiting singularities at the bulk phase 
transition in a given finite system sets in when L/£ = 0(1). 

The tacit assumption is that all other reference lengths are irrelevant and will lead only to 
corrections towards the above picture. Moreover, the crucial point in the finite-size theory is that 
we always assume: the finite-size linear L of systems under consideration and the correlation 
length £ are large in the microscopic scale. This means L ^> a and £ a, where a is the lattice 
spacing. For a comprehensive recent review on this subject and other studies related to it see Ref. 



spt 



To investigate the finite-size scaling properties of systems with long-range (LR) interaction 
decaying at large distances r as r~ d ~ a , where d is the space dimensionality and < a < 2 a 
parameter controlling the range of the interaction, one needs to extract the finite-size effects from 
the d dimensional lattice sum 



where due to the periodic boundary conditions q is a discrete vector with components qi = 2-Krii/L 
(rii = 0, ±1, ±2, •••),« = 1, • • • , d. L is the size of the box confining the system and t is a 
parameter measuring the distance to the bulk critical point. Here we will not comment on the 
restrictions on d, a and 7, nor will we dwell on the problem of convergence of (11.11) : these details 
should be clear from the context. Notice that a = 2 corresponds formally to the case of short-range 
(SR) interaction [2]. 7 is a parameter allowing to treat classical (7 = 1) and quantum (7 = 1/2) 
systems on an equal footing. Higher values of 7 appear in the investigation of finite-size systems 
to the one loop order in the field theoretical approach. For a recent review on the critical properties 
of systems with LR interaction see Ref. 0,^0. 

An other reason for considering W] a (t, L) is, as we will see below, the direct mapping be- 



tween the lattice sum (11.11) and some lattice sums in combinations with integrals that appear in the 
theory of systems with strong anisotropic LR interaction of the asymptotic form (see for example 
references lU, Id, LzL Id, IsD ) 

J(q)~ J(0)+a|[|q[||" + a ± |qj.r, (1.2) 

where the first r directions (called "parallel" and denoted by the subscript ||) are extended to 
infinity and the remaining s directions (called "transverse" and denoted by _L) are kept finite, 
with r + s = d and a±, an are metric factors and p, a > 0. Let us note that there is a limited 
number of papers that consider FSS assumptions on a microscopic models |Q,Q]. It seems that 
the considerations are mainly on phenomenological level or via computer simulations (see, for 
example refs. JvL [sL 13]) because of the problems emerging in analytical treatments. 

The study of the difference between the d-dimensional sum (11.11) at large sizes L and its limiting 
integrals is crucial in the derivation of finite-size effects. In the particular cases 7 = 1 or 1/2 to 
solve this problem several approaches have been proposed |Q, [k], [ill Q Q, Q 15. 1^. 
Am ong them the most universal one is that based on the Poisson summation formula 15, 111. 

The aim of this approach Jl^l is to factorize the (/-dimensional sum in the r.h.s 
of equation dl.lt and to reduce it to an one-dimensional effective problem. The term |q| CT in 
conjunction with 7 to be arbitrary in the interval < 7 < 00 causes peculiar mathematical 
problems concerning the evaluation of the lattice sums over q. The aim of the present study is to 
generalize the previous investigations for arbitrary < 7 < 00. By virtue of the relation 

W]+ n (t,L) = - -A—^ "wJJt,L), (1.3) 

d > a v ' ; 7(7 + 1) ■■■(j + n-l)d n t d ' aK J 

one can see that it is formally sufficient to consider only the case < 7 < 1. Let us first consider 
separately the classical and the quantum case. 
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12, 



A. Classical case (7 = 1) 



In the case of classical systems with SR interaction, corresponding to 7 = 1 and a = 2, the 
following substitution is used as an indispensable ingredient for the FSS calculations (see, e.g. 



WL(t,L) 



dx exp(— tx) 



^exp(-g 2 s) 



(1.4) 



where q is one-dimensional discrete vector. 
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This is the so called Schwinger parametric representation. The analytic properties of the func- 
tion ^ g exp(— q 2 x) are very well known, since it is nothing but the reduced Jacobi 9 3 function. 
The aim of the above procedure is two-fold: (i) to exponentiate the summand and to reduce the 
d-dimensional sum to a one-dimensional sum with well known analytic properties, and (ii) to give 
the dimensionality d the status of a continuous variable. 

In the presence of q CT term, it is not so easy to realize this procedure. The problem has been 
solved by suggesting different generalizations 11711 of the Schwinger representation 

(11.41) that lead to different obstacles. 

In order to preserve the possibility for further analytical consideration based on the properties 
of the reduced Jacobi 6 3 function in Ref. [ 1 1 ] the following representation has been used 

d 



E 



exp(— q x) 



(1.5) 



The price one pays for this is that instead of the simple exponent in the integrand of (11.41) . the 
function Q a (t) appears: 

poo 

Q a (x) = / dyexp(-xy)Q a (y), (1.6a) 



where 



QM = - 



sin 



< a < 2. 



(1.6b) 



7r 1 + 2yf cos (§7r) + y a ' 
First the connection between Q„{x) and the Mittag-Leffler type functions in the theory of FSS 
was established in reference [12]. This reads 



Q a (x) = £2 ^(-r), 
where E a ^(z) are entire functions of the Mittag-Leffler type defined by the power series 



(1.7) 



HQ 



OO I, 



k=0 



T(ak + (3) 



a,peC, Re(a) > 0. 



This shows that the study of the finite-size behaviour of the lattice sum 



Wl a (t,L)= / dxx^E.ji-txS 
Jo 



E 



exp(— q x) 



are a direct consequence of the analytical properties of E a ${z i 



Q. 



(1.8) 



(1.9) 
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B. Pure quantum case (7 = 1/2) 



For the investigation of the FSS at zero temperature when the phase transition is driven by a 
quantum parameter, in d 1 . 1 b we have 7 = |. Then the following integral representation is obtained 

1 9 r°° 1 

The auxiliary variable p 2 adds an effective extra dimension. Indeed the pure quantum system 
corresponds to a d + 1 dimensional anisotropic classical system with the geometry of a cylinder 
L d x 00. Recall that the bulk critical behaviour (e.g. critical exponents) of a pure quantum system 
is equivalent to a d + z (z = ~ dynamic critical exponent) dimensional classical system (this is the 
so called classical to quantum crossover). 

On the other hand, in the spirit of relation dl.91) the following modification has been proposed 



1 r°° 
Wl,(t,L)= / dxx^G^i-txi] 
Jo 



^exp(-g 2 x) 



The new functions G„^(z) are defined by the power series 



(lid 



G a ^z) = — V 1 J' Z - a, PeC, Re(«)>0. (1.12) 

r-i 

Some results on the analytic behaviour of these functions are presented in reference |y|]. In the 
particular case a = |, (5 — | the following identity |Q] 

G |lf (-z) = - / % f {-(z + p 2 ))dp, (1.13) 

can be obtained from Eq. (11.101) and the relation of its l.h.s and r.h.s. with the functions G^s. (z) 
and E* t z (z), respectively. 



C. Anisotropic case (0 < 7 < 1) 

In this case, instead of (11.101) we propose the following identity that can be obtained (see Ap- 
pendix |Aj) after some algebra (0 < 7 < 1) 

1 r°° 1 

(1 - 7)r(7)r(l - 7) Jo ^ t + p — + | q | ff 



5 



Equation (11.141) generalizes the result d 1 . 1 Oft corresponding to the pure quantum case. Here the 

i 

auxiliary variable p l -~i acts effectively as an anisotropic extra dimension that generates additional 
mathematical difficulties. In the denominator of the summand in the r.h.s of Eq. (11.141) one can 
easily recognize the form of the anisotropic interaction (11.21) with s = d,r = 1 and p = 1/(1 — 7). 

In this paper, we present new representation formulas for the lattice sums defined in equation 
(11.11) relevant to the investigations of the finite-size scaling properties of a large class of systems: 
classical, quantum and systems with strong anisotropy. Following the lines of consideration men- 
tioned for the particular cases of the previous subsections II Al and II Bl our aim here is to present 
functions depending on three parameters a, (3 and 7 that can play the same role as the functions 
Ea t p{z) and G a ^{z). The mathematical properties of these functions will be discussed in the next 
section and some applications will be given. 



H. GENERALIZED MITTAG-LEFFLER FUNCTIONS 

The following generalization of the Mittag-Leffler functions is defined by the power series [ 2 

are of a significant interest. Here 

(7)0 = 1, (j) k = 7(7 + 1) ( 7 + 2) ■ . - (7 + k - 1) = r ) ' k = l,2,.--. (2.2) 

These functions are named after Mittag-Leffler who first introduced the particular case with (3 = 
7=1. Recently the interest in this type of functions has grown up by their applications in some 
evolution problems and by their various generalizations appearing in the solution of differential 
and integral equations. For some mathematical applications see Refs. QQ and references 
therein. 

In the present study we will show that these functions can play an intrinsic role in the theory 
of FSS. Remark that the generalized functions (12.11) reduce to the G a> p given by (11.121) in the 
particular case 7 = ~. 

One of the most striking properties of these functions is that they obey the following identity 

POO 

(l + z)-?= / dxe~ x x l3 - 1 El l3 (-x a z), Re(7),Re(/3) > 0, \z\ < 1, (2.3) 
Jo 
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which is obtained by means of term-by-term integration of the series (I2.ll) . As we will show 
the identity (12.31) lies in the basis of the mathematical investigation of FSS in systems with LR 
interaction. If we set in the identity ( 12.31) z = y~ a , y > 0, and x = ty, we will obtain the Laplace 
transform 

„,cey— poo 

dte- yt t l3 - 1 El )/3 (-t a ) (2.4) 



from which we derive a new identity by setting (3 = cry: 

1 



dte -vt t <rr-iE7 (_*«). (2.5) 



(i + y a V Jo 

With the help of the above identity one immediately obtains the relation 

d 

Re( 7 ) > 0, (2.6) 



poo 

f a (t,L) = / dxx^El J-tx^ 
Jo 2 1 



^exp(-g 2 x) 



/oo 
Elz (-{z + p^)j dp, < 7 < 1. (2.7) 



which is the quested generalization of ( 11.91 ) and (11.1 II) . Now it is easy to obtain from Eqs. ( 11.91 ). 
(11.141) and (12.61) the generalization of the integral relation dl.131) : 

^ ( "" ) = (i-7)r(7)r(i- 7 ) Jo 

Notice that El z := £?|,f • The asymptotic expansion for z ^> 1 of the generalized Mittag-Leffler 
functions E^ Jz) can be obtained from (see Appendix iBl): 

In the particular case /3 = 0:7, relevant to the physical situations we are discussing here, equa- 
tion (12.81) reduces to 



El, ai {-x) ~ -7-f^-v, x » 1, (2.9) 



r(-a) 

where only the leading term is accounted for. 

In the following section we will discuss some applications and relations to previous results 
obtained in the framework of the FSS investigations in the classical and the quantum cases. 

IE. FINITE-SIZE COMPUTATIONS 

For FSS computations of 0(n) systems one needs the large-L behaviour of normalized lattice 
sums defined by 

r ^ i »4^FW' (31) 

q^O V ^ lHl I 



where 7 is relevant to different physical cases. Indeed for integer 7 we have classical systems, 
while for half-integers we have the quantum situation and for 7 < 1 systems with geometry 
L r x 00 s and strong anisotropy of the type (I1.2I ). The last follows directly from the identity 
(compare with (11.141 )') 



1 



r(i_i)r(i; 



t + pf + Iql' 7 



p>\. 



(3.2) 



The method we use here to extract the large-L behaviour of (13.11) is based upon the identity 
(12.6b . After rearrangement of d3.ll) we obtain 



jjya—d poo 



(2< 



7<T 



dxx 1 " 2 1 E 1 



where 



A(x) = e ~ 



tL° 



(2vr) CT 



X 2 



[A d {x) - 1] , (3.3a) 



(3.3b) 



For large x, A(x) — 1 decreases exponentially and the integral in the right-hand side of equation 



(I3.3al) converges at infinity. For x — ► 0, the Poisson transformation formula 

r2~ 



A(x) 



A 



(3.4) 



shows that A(x) converges. 



For small x the integral in the right-hand side of equation (I3.3al) has an ultraviolet divergence 
for Re(d) > 70-. So, an analytic continuation in d is required to give a meaning to the integral. 
Adding and subtracting the small behaviour of the function A(x), we get after straightforward 
algebra, 

m a (t, L) = L~ d+ ^ \m {tL*)^ + F] (tV) 



where the constant 



and the functions 



n 1 



2 1 r( 7 -j)r(g) 
a( 47r )f r( 7 )r(f) 



(2tt)^ Jo 



dxx 1 2 E 



1 p7 



(2vr) CT 



X 2 



A d (x) 



7T \ 2 

X 



(3.5a) 



(3.5b) 



(3.5c) 



The first term in (I3.5al) is the bulk contribution (it is L-independent) and the second term is the 
corresponding finite-size correction. The form (I3.5al) is suitable for the investigation of FSS in 
the vicinity of the critical point i.e. t ~ 0. The function Fj a (y) enters in the expressions for the 
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scaling functions of various thermodynamic observabales. The dependence on the linear size L of 
the system is included in the scaling variable y. The behaviour of any thermodynamic function 
is tightly related to the asymptotic behaviour of Fj a (y), which in turn depend upon that of the 
Mittag-Leffler functions. For detailed discussions of different models with the particular values 
7 = 1 and 7 = \, the reader is invited to consult references , where the finite-size 

scaling predictions are investigated in great details. Let us not that at this level, the anisotropy of 
the scaling behaviour in Eq. d3.5al) only appears through the parameter 7. 

By setting t = in (I3.5al> we obtain an expression for the finite-size shift of the bulk critical 
parameter driving the phase transition. This is proportional to Fj (T (0)L~ x , where A = d — 70" 
is the shift critical exponent for the specific value of 7. The coefficient Fj a (0) can be evaluated 
analytically as well as numerically for different values of the free parameters d, a and 7 using the 
method developed in reference [25]. 

According to the standard FSS we must have A = l/u, where \ jv is the critical exponent 
measuring the divergence of the correlation length. The value of v depends on the concrete micro- 
scopic model. For illustration, in the particular case of symmetric 0(n) model in the limit n — > 00, 
one can consider two cases: classical and quantum, where v — l/(d — a) and v = l/(d — cr/2), 
respectively. In both cases our result confirms the FSS theory predictions (see,e.g. |Q]). Fur- 
thermore in order to make contact with the case of strong anisotropy of refs. 1^,13] the effective 
dimensionality D = 2d /a + 2/ p must be introduced. It determines the conditions for the phase 
transition to take place, if the anisotropic LR has the asymptotic form (11.21) with s — 1, r — d. In 
this case v = 2/(a(D — 2)) and again we have agreement with FSS theory , provided 2 < D < 4. 
Moreover, introducing the effective dimension .0(7) = 2d/a + 2(1 — 7) we can consider a more 
general classical system that includes 7 = 1, 1/2, (1 — 1/p) as particular cases. 

The finite-size correction to the bulk critical behaviour ( of e.g. susceptibility) can be extracted 
from the asymptotic form of the functions Fj a (y) defined by (I3.5cl) at large argument. This in turn 
can be obtained with the help of the expansion (12.9b . After some algebra we get (see Appendix ICl) 



for a = 2 i.e. for the SR case and 



FJM 



-y 



2^(27r)W( 7 ) 



i(d-2 7 -l) -y/y 



(3.6a) 



v 2) 1^0 11 



-(1+7) 



(3.6b) 



9 



for < a < 2, corresponding to the LR case. Equations (13.61) generalize equations (3.29) of 
reference IieU obtained for the particular case j — 1. 

Equations (13.61 ) reflect the fact that for systems with LR interaction the finite-size corrections 
fall-off in power law rather than exponential as it is the case for their counterparts with SR interac- 
tion. These results are generalizations of those obtained previously in the case of classical systems 



and those obtained for their quantum counterparts Q, 



1511 . The former cases can 



be obtained by using integer values for 7 and the latter ones by using half-integer values. 
IV. CONCLUSION 

We presented some mathematical results on the investigation of the FSS in 0(n) systems based 
on the generalized Mittag-Leffler functions (12.11) that have well known analytic properties. Mainly 
two type of systems are of particular interest. 

(i) The fully finite (/-dimensional systems with LR interaction decaying algebraically with the 
interparticle distance. 

This is the case with < a < 2 in (ll.lt . A special emphasis on the mathematical difficulties 
arising in the investigation of the FSS both in the classical (with 7 = 1 in (11.11) ) and quantum 
cases (with 7 = 1/2 in dl.lt ) are discussed. The used techniques allow the investigations to be 
simplified and express the results for various thermodynamic quantities in terms of simple, with 
well defined analytic properties, mathematical functions. An integral representation (12.61) to deal 
with such difficulties, at least asymptotically, are presented. It turned out that both cases can be 
treated on an equal footing. 

(ii) The classical system with mixed geometries with both finite and infinite sizes and strongly 
anisotropic interaction of the type dl .2ft . 

Such type of systems are considered in reference [6], where < p, a < 2 and 7 = 1 — 1/p. An 
other interesting case is the m-fold Lifshitz point that is characterized by an instability associated 
with the absence of quadratic terms in the form q 2 a in the effective Landau-Ginzburg-Wilson 
Hamiltonian for all a = 1, 2, • • • , n < d O, 3- Then in (II .2I> one must set p = 4, a = 2 and 
7 = 3/4. This gives a simpler way of solving such problems using generalized Mittag-Leffler 
functions. It is based on the established mapping d3.2ft to a fully finite-size system with specific 
< 7 < 1 in (fO . 

In conclusion, our considerations establish that we can study finite-size scaling behaviour of 
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classical systems, quantum systems and systems with strong anisotropy confined in mixed geome- 
try oo x L d , in the framework of a fully finite anisotropic system with a classical critical behaviour. 
This is achieved in a unified fashion, varying the superscript 7 in the generalized Mittag-Leffier 
functions. 
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APPENDIX A: DERIVATION OF EQUATION (fTTHA 

We have the relation 



lp-«/i* T f-\= [ x a - l e- pxtl dx- n, Re(a), Re(p) > 0. (Al) 
V V/V Jo 

or if = 7 

= r( 7 )7 



= ^ / x- {1 ~^e~ px "dx. (A2) 



Using twice (IA2I) we get 

1 [A v 



p-r r( 7 )r(i- 7/ i)7 J 

Now on the free parameters \x and v we impose the conditions 



roc poo 

dx dte'^e-^t^-^- 1 . (A3) 
Jo Jo 



11-7^-1 = 0; A* = 1 (A4) 



and obtain the identity (7 < 1) 

1 1 r° . 1 



P 7 (1 _ 7 )r( 7 )r(i - 7 ) y p + t jh 

Eq. (11.141) immediately follows from the above identity. 



dt — . (A5) 



APPENDIX B: DERIVATION OF THE ASYMPTOTIC BEHAVIOUR OF THE GENERALIZED 
MITTAG-LEFFLER FUNCTIONS <2jJ 

An integral representation of the generalized Mittag-Leffier functions E^Jz) can be obtained 
with the aid of the Henkel integral for the inverse gamma function 

1 1 



T(z) 2m 
11 



e u u~ z dz, (Bl) 



c 



where the integration contour C is a loop which starts and ends at x = — oo and encircles the 
origin in the positive sense: — n < arg z < n on C [12811 . This enables to get the result 

1 



2ma 



dv 



e vV a v -y-l+(l-P)/a 



(B2) 



c 



(v - zp 

In the following we will investigate the asymptotic behaviour of the generalized Mittag-Leffier 
functions at large argument, following the method used in reference [5]. This may be performed 
with the aid of the series [29] 



(x + z) 7 = X 7 



fc=0 



k\ 



x 



r. 


z 




X 



z 

X 



(B3) 



After substitution of the latter equation into the integral representation (IB2I) one obtains (12.81) . 
APPENDIX C: LARGE ASYMPTOTIC BEHAVIOUR OF F] a (y) FROM @3cJ 



To obtain the large y asymptotic behaviour (13.61) of the functions F] a (y) we rewrite (I3.5cl) . 
with the help of the identity (13.4ft . in the form 



d 
712 



dxx 1 ^ 2 1 E~'l c_ 

2 ' ' 2 



-y 



(2vr) 7CT J 



dxx 1 ' 2 l EZ a I — y 



X2 

(2tt) ct 

a_ 

xi 



A d — 



(CI) 



Using the identity 



(C2) 



/•OO 

/ dxx^EZ „ (-xt) =1, a > 
Jo 2 ' 2 

from the second term of equation ( IC1I ) we obtain the first terms of equations (I3.6I ). 

Further, taking into account the asymptotic behaviour (12.81) of the functions E^Jz) and after 
subsequent integration in the first term of equation (IC 1 ft . we obtain finally the asymptotic behaviour 
given by equations (13.6ft . 
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